Nontrivial Exponent for Simple Diffusion 
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The diffusion équation dt<f> = X7 2 (f> is considered, with initial condition 4>(x, 0) a gaussian random 
variable with zéro mean. Using a simple approximate theory we show that the probability p n (pi,t2) 
that </>(x, i) [for a given space point x] changes sign n times between t\ and tu has the asymptotic 
form p„(t 1 ,t 2 ) ~ [In(t2/ti)] n (ti/t2r 9 . The exponent 6» has predicted values 0.1203, 0.1862, 0.2358 
in dimensions d = 1, 2, 3, in remarkably good agreement with simulation results. 



The diffusion équation, d t (f> = V 2 </>, is one of the 
fundamental équations of classical physics. The exact 
solution of this simple équation, for an arbitrary ini- 
tial condition </>(x, 0), can be written down explicitly: 
<j>(x,t) = Jd d xG(x ~ x',i)0(x',O), where G(x,i) = 
(47rf)~ d / 2 exp(— x 2 /4t) is the Green's function (or 'heat 
kernel') in d dimensions. The solution is characterized 
by a single growing length scale, the 'diffusion length' 
L(t) ~ i 1 / 2 . It may corne as a surprise, therefore, to dis- 
cover that there is a nontrivial exponent associated with 
this simple process. 

It is the purpose of this Letter to point out that the so- 
lutions of the diffusion équation exhibit some remarkable 
and unexpected properties associated with their time 
évolution, and to présent a simple theory which accounts 
for this behaviour. We consider specifically a class of ini- 
tial conditions where </>(x, 0) is a gaussian random vari- 
able with zéro mean. Our basic question is the following. 
What is the probability po(t) that the field <f> at a partic- 
ular point x has not changed sign up to time i? Précise 
numerical simulations in d = 1 and 2, discussed below, 
demonstrate a power-law decay of the form po(t) ~ t~ , 
with = 0.1207±0.0005ford = 1, and 0.1875±0.0010for 
d = 2. We will présent a simple analytic treatment which 
gives results in extraordinarily good agreement with the 
simulations. Furthermore, the analysis gives the more 
gênerai resuit p n (ti,h) ~ [ln^Ai)]™ {ti/t 2 )~ e for the 
probability that the field changes sign n times between 
t\ and t2, for t 2 3> t\. The key idea underlying thèse re- 
sults is that the gaussian process </>(x, t) is a gaussian sta- 
tionary process in terms of a new time variable T = hit. 
The central assumption in the analysis is that the inter- 
vais between successive zéros of </>(x, T) can be treated 
as independcnt. 

Exponents analogous to that introduced above havc 
recently excited much interest in other contexts [p]dlÔ|. 
The simplest such System is the d = 1 Ising model at 
température T = 0. For évolution under Glauber dy- 
namics from a random initial state, the probability that 
a given spin has not flipped up to time t decays as t~ e , 
with = 3/8, though the proof of this is surprisingly sub- 
tle H . This d = 1 method is difficult to extend to higher 
dimensions, although values for have been obtained nu- 
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dimensions has recently been developed [M , whose prédic- 
tions are consistent with simulation results. In gênerai, 
the non-triviality of po (t) is a conséquence of the fact that 
it probes the entire history of a non-Markovian process. 

We begin by presenting the theoretical approach and 
the numerical simulation results. Expérimental ramifica- 
tions will be discussed briefly at the end of the Letter. 
Other contexts in which a nontrivial exponent might 
be expected will also be discussed. 

The starting point for the discussion of the diffusion 
équation is the expression for the autocorrélation func- 
tion of the variable X(t) = 0(x, t)/([</>(x, t)] 2 ) 1 ' 2 for 
some fixed point x. For 'white noise' initial conditions, 

x' 



(0(x,O)0(x',O)) -5 d (x 

a{h,t 2 ) = (X(ti)X(t 2 )) = [4tita/(ti 



this takes the form 
t 2 ) 2 ] d '\ 



(1) 



More generally, this form is asymptotically correct pro- 
vided the initial condition correlator is sufHciently short- 
ranged (it must decrease f aster than |x — x'|~ d ). 

Introducing the new time variable T = hit, one sees 
that the autocorrélation function becomes a(Ti,T2) — 
/(Ti - T 2 ), where f(T) = [sech(T/2)] d / 2 . Thus the 
process X(T) is stationary (the gaussian nature of the 
process ensures that ail higher-order correlators are also 
time-translation invariant). This is an important simpli- 
fication. Note that the anticipated form of the probabil- 
ity of X(t) having no zéros between ti and t 2l Po(ti, t 2 ) ~ 
{t\/t 2 ) for t 2 t\, becomes an exponential decay, 
Po ~ exp[— 0(T 2 — Ti)], in the new time variable. This 
reduces the calculation of an exponent to the calculation 
of a decay rate (7). The only approximation we shall 
make is that the intervais between successive zéros of 
X(T) are statistically independent. This 'independent 
interval approximation' (IIA) was introduced in another 
context some forty years ago |ÏÏJ . We shall find that it is 
an extraordinarily good approximation for the diffusion 
équation. 

As a preliminary step, we introduce the 'clipped' vari- 
able a = sign(X), which changes sign at the zéros of 
X{t). Clearly, the correlator A(T) = (ct(O)ct(T)) is de- 
termined solely by the distribution P(T) of the intervais 
between zéros. The strategy is to détermine P(T) from 
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A(T), and p (T) from P(T). To this end we note first 
that 



A(T) 



(a(T)] = hm' 1 {[sech(T/2)] d ^), 



(2) 



where the first equality holds for any gaussian process. 

Next one expresses A(T) in terms of the interval-size 
distribution P(T). Clearly 



A (T) = V(-1)>„(T), 



n=0 



(3) 



where p n (T) is the probability that the interval T con- 
tains n zéros of X(T). We define Q(T) to be the prob- 
ability that an interval of size T to the right or left of 
a zéro contains no further zéros. Then P(T) — —Q'(T). 
For n > 1 one obtains immediately 

Pn (T) = (T)- 1 f dT x [ dT 2 ... f dT n x (4) 

xQ(T x )P{T 2 - Ti) . . . P(T„ - T n _ x )Q{T - T n ), 

where (T) is the mean interval size. One has made 
the IIA by writing the joint distribution of n succes- 
sive zero-crossing intervais as the product of the distribu- 
tion of single intervais. Taking Laplace transforms gives 
Pn(a) = [Q(s)] 2 [P( S J]"-V<r). But P(T) = -Q'(T) im- 
plies P(s) = 1 — sQ(s), where we have used Q(0) = 1. 
Using this to eliminate Q(s) gives the final resuit 



Pn(s) 



'l-P(s)) (P(s)) , n>l, (5) 



(T)s 2 



(G) 



where the resuit for p~o(s) follows from the normalization 
condition X^ p n (*) = ^ wnic h g ives S,T=iP"( s ) = 

1/S - 

Finally the Laplace transform of (Q) gives A(s) = 
S^Lo( — l)™Pn( s )- Performing the sum employing (||) 
and (||), and using the resuit to express P(s) in terms 
of A(s) gives the desired resuit 



where 



P(s) = [2-F( S )}/F(s), 



F(s) = l + ((T)/2) S [l- S Â(s)]. 



(7) 



(8) 



Equations (|||]) are a gênerai conséquence of the in- 
dependent interval approximation. The function F{s), 
defined by (||) , is completely determined by the autocor- 
rélation function A(T), and contains ail the information 
needed to compute the probabilities p n (T). We have in 
mind, of course, to apply this approach to the diffusion 



équation, where A(T) is given by (||). For this case the 
mean interval size (T), required in (||), can be simply 
evaluated. For T — > 0, the probability to find a zéro 
in the interval T is just T/(T), so A(T) ->• 1 - 2T/(T). 
This gives (T) = -2/A'(0) = ir^/S/d, using @ in the 
final step. 

We note a very important point at this stage. The fact 
that A'(0) is finite (i.e. /'(0) = and /"(0) ^ 0) is spé- 
cial to the diffusion équation, which allows us to use the 
IIA. Physically this means that the density of zéros is a 
finite number. However, for many Gaussian stationary 
processes, such as the one that arises in an approximate 
treatment of the Ising model 0, /'(0) ^ 0, implying that 
A'(0) diverges. In this case, the IIA cannot be used. For 
such processes, the variational and perturbative methods 
developed in Réf. give reasonably accurate results. 

The asymptotics of P o(T) are controlled by the singu- 
larity of po{s) with the largest real part, i.e. [from (0)] by 
the corresponding singularity of P(s). The expectation 
that po(T) ~ exp(— 9T) suggests that this singularity is 
a simple pôle, i.e. that F (s) has a simple zéro at s = —9. 
Using (g) in (m), and inserting (T) = ny/8/d, gives 



F (s) = !+,-[- 



1/2 



2 S f°° 
1 / dT exp(-sT) x 

T JO 



x sm 



srch"' 1 { | 



(9) 



Clearly F(Q) = 1, while F(s) diverges to — oo for s — > 
—d/A. Between thèse two points F (s) is monotonie, im- 
plying a single zéro in the interval {—d/A, 0). Solving (||) 
numerically for this zéro, and identifying the resuit with 
—9, gives the values of 9 shown in table 1. For future 
référence, we note from (0) that the residue R of the cor- 
responding pôle of P(s) is R — 2/F'(— 9). The values of 
R, which controls the amplitude of the asymptotic de- 
cay of p n (T), are also given in table 1. Recall that the 
behaviour po(T) ~ exp(— 9T) translates in 'real' time to 
a decay law po(t) ~ t~ e for the probability that cf> at a 
given point has not changed sign. It is also easy to ex- 
tract the large-ci behaviour of 9 from Eq. (9): we find, 
to leading order in d, 9 sa 0.145486Vd. 

The predicted values of 9 were tested in d = 1 and 
2 by numerical simulations. The diffusion équation was 
diserctized in space and time in the form 



i(t + 1) = 4>i{t) + a^ifait) - 



(10) 



where j runs over the nearest neighbours of i on a linear 
(d = 1) or square (d = 2) lattice. A stability analysis 
shows that the solution is unstable for a > a c — l/(2d). 
Preliminary studies showed that the asymptotic expo- 
nent is independent of a for a < a c , but that a value 
a = a c /2 seems to give the quickest onset of the asymp- 
totic behaviour. This value was therefore used in ail sim- 
ulations reported here. Systems of 2 20 (2 22 ) sites were 
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studied in d — 1 (d — 2), for times up to 2 17 (2 12 ). Data 
for longer times in d = 2 suffer from noticeable finite-size 
effects. The initial values of (\>i were chosen indepen- 
dently from a gaussian distribution of zéro mean. Using 
a rectangular distribution gave the same asymptotic ex- 
ponent within the errors. Several random number gener- 
ators were tried: Ail gave consistent results (within the 
errors). 

The simulation results are presented in Figure 1. The 
data are an average of 17 [d — 1) and 22 (d = 2) runs 
with independent initial conditions. An effective expo- 
nent 9{t) is extracted from a least-squares fit of log 2 po 
against log 2 t over 5 consécutive values of log 2 t. The 
error bars shown in the figure were obtained from the 
fits. The resulting exponent 6(t) is then plotted against 
1/ log 2 t, where here log 2 t is the mid-point of the 5 val- 
ues. The best estimâtes of 9, shown in table 1, were 
obtained by plotting t 6 po (t) against log 2 t and choosing 
9 such that, after an initial transient, the data show no 
systematic upward or downward trend with increasing t. 
The agreement with the theoretical prédictions (table 1) 
is quite remarkable, showing that the IIA is an extraor- 
dinarily good approximation in this context. 

The case of correlated initial conditions is also of inter- 
est. If the Fourier-space corrélations are (0k(O)0-k(O)) ~ 
k a for k — > (a > —d), the autocorrélation function of 
X(t) still has the form (|ï|), but with d replaced by d + a. 
The independent interval approximation, therefore, pre- 
dicts that the dependence of 9 on d and a enters only 
through the combination d + a. To test this, we sim- 
ulated the case d = 1, a = 2, noting that a = 2 corre- 
sponds in real space to differentiating uncorrelated initial 
conditions (or taking finite différences on a lattice). The 
resuit from 12 runs (Figure 1) is 9 = 0.2380 ± 0.0015, 
close (as anticipated) to the predicted resuit 0.2358 for 
uncorrelated initial conditions in d = 3. 

The asymptotics of the probability p n (t\, t 2 ) for having 
n zéros between times t\ and t 2 are also readily calcula- 
ble within the IIA. From (H|) and (|J), the singularity in 
Pn(s) as s — —9 is an (n + l) th -order pôle of strength 
R n+1 1 (T)9 2 , where R is the strength of the simple pôle 
in P{s). Inverting the Laplace transform, and retaining 
only the leading large-T behaviour, gives (for ail n) 



Pn(T) 



R (RT) r 



(T)9 2 



exp(-é»T). 



(H) 



With T = \n(t 2 /ti), one obtains 

Pn{hM) - (R n+1 /(T)9 2 ) [Ht 2 /h)] n (h/t 2 ) e . (12) 



finds m p n {t) ~ (n)-V(n/(n)), with (n) ~ y/t. [The 
exponent 9 in those Systems émerges from a singular be- 
haviour of the scaling function f(x) as x — > 0. Note that 
in the présent work (n) ~ T ~ Int.] 

We turn to a brief discussion of the expérimental rele- 
vance of our results. The ubiquity of the diffusion équa- 
tion in physics implies that applications will be many 
and varied. As a concrète example, however, considcr 
the reaction-diffusion process A + B — * C, where C 
is inert and immobile. The corresponding rate équa- 
tions for the concentrations are dnj^jdt = ~V 2 ua — R, 
dus/dt — V 2 hb — R, and dnc/dt — R, where R is the 
reaction rate per unit volume (R oc uatib for d > 2 p3|). 
The concentration différence, An = ua - îIb, obeys the 
simple diffusion équation. If the A and B species are ran- 
domly mixed at t = the System evolves, for d < d c = 4, 
to a coarsening state in which the two species segregate 
into domains |Ï4| , separated by domain walls whose loca- 
tions are defined by An = 0. Subséquent production of 
the inert species C is slaved to the motion of the domain 
walls, which are zéros of the diffusion field An. The frac- 
tion of space not infected by the C species will therefore 
decay asymptotically as t~ e . 
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FIG. 1. Effective exponents 8(t) plotted against l/log 2 t 
for the diffusion équation in d = 1 (lower data set), d — 2 
(middle set), and d = 1 with correlated initial conditions (up- 
per set). The downturn in the upper set at late times is not 
statistically significant (note the larger errors on the last two 
points). The best estimâtes of are given in table 1. 



When the time t\ corresponds to the initial condition, one 
has to set ti equal to a constant of order unity, as was im- 
plicit in the earlier treatment of Poit). Setting t% = t one 
then gets p n (t) ~ (lnt)™ t . This rather strange-looking 
resuit does not have the scaling form found in the voter 
model and in Ising Systems in d — 1 and 2, where one 



We conclude with other examples of non-trivial expo- 
nents 9 which have not been addressed in the literature. 
The first is associated with the dynamics of the global 
order parameter M(t) (e.g. the total magnetization of 
an Ising ferromagnet) at a critical point T c , following 
a qucnch to T c from the high-temperature phase. The 
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quench prépares the System in a state with random initial 
conditions. In the subséquent évolution (now stochastic, 
rather than deterministic), the probability that M if) has 
not changed sign since t = decays as t~ 9c , where 9 C is 
a new critical exponent |Ï5[ . For reasons similar to those 
given for the diffusion problem, we expect 9 C to be an in- 
dependent exponent, i.e. not related by any scaling law to 
the usual static and dynamic exponents. As a second ex- 
ample, one can consider M(t) for a quench to T = from 
high température. In this case, po(i) ~ t~ 00 , where Oq 
differs from the corresponding exponent for single spins. 
For the d = 1 Glauber model, for example, the prob- 
ability that M(t) has not changed sign decays with an 
exponent 9q = 1/4 which differs from the exponent 
3/8 obtained for the zero-flip probability of a given spin 



As a final example, consider the generalised one- 
dimensional random- walk équation d n x/dt n = £(t), 
where £ is gaussian white noise. The cases n = 1,2,... 
correspond to a random velocity (the usual random 
walk), random accélération, etc. The first two 9 n are 

91 = 1/2 and 9 2 — 1/4 [0, but larger n have not 
been considered before to our knowledge. Application of 
the independent interval approximation |Ï7j| gives équa- 
tions of the same structure as for the diffusion pro- 
cess, but with sech d (T/2) in @ and (||) replaced by 
(2 - l/n) exp(-T/2) 2 F X [1, 1 - n; 1 + n; exp(-T)], where 
2-F\ is the hypergeometric function. This approach gives 

9 2 = 0.2647 (instead of 1/4) while, for larger n, 9 n ap- 
proaches a limiting value 9^ — 0.1862 . . ., i.e. the same 
exponent as the d = 2 diffusion équation! In fact, the 
equality of the exponents for the n — oo process and 
d = 2 diffusion can be proved exactly $tï\ , implying a 
limiting exponent 0.1875 ± 0.0010 (from table 1) for the 
former. 

To summarize, we have calculated the probability for n 
zéro crossings, between times ii and < 2 , of a diffusion field 
at a given point in space, by assuming that the intervais 
between crossings, measured in the variable T = hit, are 
independent. The time-dependence of thèse probabili- 
ties is characterized by a single non-trivial exponent 9, 
the predicted values of which are in excellent agreement 
with précise simulation results in one and two dimen- 
sions. Thèse ideas are relevant to any System where the 
diffusion équation (or 'heat équation') plays a rôle, rang- 
ing from physical and chemical Systems to fluctuations in 
financial markets, and can be extended to other gaussian 
processes. 
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0.1203 
0.1862 
0.2358 
0.2769 
0.3128 



0.1207 ± 0.0005 
0.1875 ± 0.0010 
0.2380 ± 0.0015* 



0.1277 
0.2226 
0.2940 
0.3527 
0.4033 



TABLE I. Exponents from theory (6th) and simulations 
(Osim), and the value of the residue R (see text), for various 
spatial dimensions d. The l d = 3' simulation resuit (*) refers 
to a d = 1 simulation with correlated initial conditions (see 
text). 
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